Orbital polaron in double exchange ferromagnets 
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We investigate the spectral properties of the two-orbital Hubbard model, including the pair hop- 
ping term, by means of the dynamical mean field method. This Hamiltonian describes materials in 
which ferromagnetism is realized by the double exchange mechanism, as for instance manganites, 
nickelates or diluted magnetic semiconductors. The spectral function of the unoccupied states is 
characterized by a specific equidistant three peak structure. We emphasize the importance of the 
double hopping term on the spectral properties. We show the existence of a ferromagnetic phase 
due to electron doping near n = 1 by the double exchange mechanism. A quasi- particle excitation 
at the Fermi energy is found that we attribute to what we will call an orbital polaron. We derive 
an effective spin-pseudospin Hamiltonian for the two-orbital double exchange model at n = 1 filling 
to explain the existence and dynamics of this quasi-particle. 

PACS numbers: 71.27. +a, 71.10.Fd, 75.25. Dk 



I. INTRODUCTION 

The double exchange mechanism for ferromagnetism 
in magnetic perovskite compounds was proposed in 1951 
by Zener^i and further developed by de Gennes^ It is 
a mechanism that works for transition metal ions with a 
partially filled d-shell building up a local spin. Additional 
charge carriers in the d-orbitals that move through the 
crystal prefer a ferromagnetic arrangement with the local 
spin according to Hund's first rule and induce in this 
way a ferromagnetic long range order. This mechanism is 
the origin of ferromagnetic behavior in manganites which 
became of technological importance after the discovery 
of giant magnetoresistance (GMR) in these compounds. 
This discovery led to a renew of interest in the double 
exchange mechanismi 3 - 

Manganites are only one class in the fascinating world 
of transition metal compounds that also include nick- 
elates, cuprate superconductors and its parent com- 
pounds, and the recently discovered iron pnictide super- 
conductors. The common feature to all these compounds 
is the strong correlation in the d-shell which leads to bad 
metal behavior, magnetism, spectral weight transfer, su- 
perconductivity, orbital excitations and other phenomena 
that are far from being well understood. To clarify the 
role of orbital fluctuations (that are absent in cuprate su- 
perconductors) we consider here the minimal model with 
two orbital degrees of freedom per site. We consider the 
region of strong electron correlations, i.e. large values of 
the on-site Coulomb repulsion U. The system is insu- 
lating at quarter filling (n = 1) and we will show that 
additional electrons (n > 1) lead to ferromagnetism by 
the double exchange mechanism that we investigate in 
detail. 

Other materials where the double exchange mechanism 
may appear are diluted magnetic semiconductors (DMS), 
a prominent example is Mn-dopcd GaAs. In this system, 



however, the kinetic p-c?-exchange mechanism, is gener- 
ally believed to be more important^ This mechanism, 
also, proposed by Zener in the same year (1951), demands 
charge carriers (usually holes) in the p-orbitals. Although 
one may expect the kinetic p-d exchange mechanism in 
GaAs:Mn, in GaN:Mn the double exchange mechanism 
is likely to be more important. The reason is the lo- 
calization of the Mn 3d orbital which is much more pro- 
nounced in GaN:Mn than in GaAs:Mni£ Despite the con- 
siderable interest and effort put into unraveling the in- 
teraction mechanisms in DMSs, no clear criteria to dis- 
tinguish between both Zener's mechanisms is known. In 
the following, we analyze the double exchange mecha- 
nism and propose a characteristic that might be used to 
distinguish it from other mechanisms. We report the dis- 
covery of a narrow quasi-particle peak at the Fermi level 
with a spin parallel to the ferromagnetic order. 

Narrow quasiparticle peaks corresponding to large ef- 
fective masses are already well known in the physics of 
strongly correlated electrons. The most prominent exam- 
ples are the Kondo peak of heavy fermions and the spin 
polaron in the t-J model for cuprate superconductors^— 
In both cases, there are spin fluctuations coupled to the 
moving charge leading to the heavy mass, and in both 
cases, the spin correlation between the local magnetic 
moment and the itinerant charge carrier is of antiferro- 
magnetic nature. We will show that both aspects are 
different for quasiparticles in double exchange ferromag- 
nets: the electrons are coupled to orbital fluctuations 
and the arrangement between local spin and itinerant 
electron is ferromagnetic. In our case, spin fluctuations 
are excluded since the system is close to saturated ferro- 
magnetism (n = 1 and \M\ = 1). But at quarter filling, 
the system shows alternating orbital order. Therefore, 
there exists a close analogy between a hole moving in 
the antiferromagnetic background of cuprate supercon- 
ductors (spin polaron)£~— and an electron moving in the 
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alternating orbital environment of double exchange fer- 
romagnets. We will use this analogy to present a concise 
interpretation of the quasi-particle that we will call an 
orbital polaron. The term "orbital polaron" was orig- 
inally introduced by R. Kilian and G. Khaliullin for a 
quasi-particle for which the charge degree of freedom is 
not only coupled to orbital fluctuations, but also to the 
latticed Later the term was used in studies on effective, 
low-energy t-J like Hamiltonians in the field of mangan- 
ites for orbital quasi-particles, but without coupling to 
the lattice^ - — We re-discover the orbital polaron here 
starting from the complete two-orbital Hubbard Hamil- 
tonian with the pair-hopping term. 

The double exchange mechanism in the two-orbital 
Hubbard model was already studied before. Many works 
concentrated on the case n = 1. Below, we use the dy- 
namical mean field theory (DMFT) with the non-crossing 
approximation (NCA) as an impurity solver. The phase 
diagram in DMFT was already calculated in Refs. [ID 
and [HI, however, the paired hopping between different 
orbitals was neglected. On the other hand, the spec- 
tral weight transfer of the complete Hamiltonian was 
analyzed in an analytical way and several characteristic 
features were found that we find back using the DMFT 
approach^ Here, we include the pair-hopping term (or 
double hopping term) into our study and we concen- 
trate on the spectral function in the slightly doped case 
(n > 1). We will show that the pair-hopping term has 
indeed only small influence on the phase diagram, but 
it is very important for the spectral function, ft leads 
to characteristic features in the unoccupied parts of the 
spectral function confirming the analytical results of Lee 
and Phillipsjii 

The paper is set up as follows. After introducing 
the double exchange Hubbard model and the numerical 
method (Sees. |H] and . we present the magnetiza- 
tion diagram in Sec. HVl Special care is taken to ana- 
lyze the spectral properties and the importance of the 
pair-hopping term is emphasized. In Scc. lIVCl thc quasi- 
particle peak at the Fermi level is investigated. To under- 
stand the model at quarter filling (n = 1) we derive an ef- 
fective low-energy model in the large U-c&se in Scc.fVland 
we show that it leads to a ferromagnetic phase with alter- 
nating orbital order. That observation allows a concise 
interpretation of the sharp quasi-particle peak in terms 
of an orbital polaron fSec. IVI|) . The paper is closed with 
a discussion and conclusion (Sec. IVII[I . 
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where r is the site index, a, /? = a, b is the orbital num- 
ber, a =y, \. is the spin projection, g is a vector that joins 
nearest neighbors. The operators S r . a = (S x , S y , S z ) rtCt 
are the spins of an electron on site r in orbital a. To find 
cigenstates of the spin term it is convenient to rewrite 

it S r - t a ■ §r,b = \ \^t,a^r,b ^r.a^^bj ^~ ^r,a^r,b' ^he 

Hamiltonian ©-(HI) describes Coulomb interaction Hd 
in Kanamori's simplified form, where the parameter U 
denotes the repulsion of electrons occupying the same or- 
bital, V = U — 2 J is the repulsion of electrons occupying 
different orbitals, J is the Hund exchange term, J' = J is 
the double hopping term, which is always present when 
exchange interactions take place. Here we retain differ- 
ent notations for the parameters J, J 1 in order to clarify 
the origin of different terms, and to facilitate the com- 
parison with the results of some papers (e.g. ReffLsh 
where J' is not taken into account. The term H t is the 
kinetic energy. For simplicity we only introduced hop- 
ping terms that are diagonal in orbital and spin indices 
and we consider a semi-circular density of states (DOS) 
being realized on the Bethe lattice. The band width is 
W = At — 2 and we chose the half-bandwidth as the unit 
of energy, which is of the order of the electron volt. We 
take the Coulomb repulsion energy U = 4W = 8 and the 
Hund exchange term 0.5 < J < 2, which is a relatively 
small parameter. 

To solve this Hamiltonian numerically we use dynam- 
ical mean field theory with the non-crossing approxima- 
tion. The following paragraph gives a quick introduction 
to this method. 



III. NUMERICAL METHOD 



II. DOUBLE EXCHANGE HUBBARD MODEL 

We will write the double exchange Hamiltonian as H = 
H t + H d , with 

H t = t ^ S rol (; d r j rg>a ^, (1) 



In the framework of the dynamical mean field theory 
(DMFT)i2, all fcrmionic degrees of freedom except those 
for a central site r = are integrated out, and the double 
exchange Hamiltonian described above can be mapped 
onto an effective single impurity Anderson model (SIAM) 
described by an effective Hamiltonian H e : 

H cS = Hf + Hf , (5) 
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where 
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refers to the impurity, which is coupled to the effective 
medium 
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represents the hybridization between the site r = 
and the effective medium corresponding to orbital a. 
is the band energy of the corresponding effective 

(respectively c a % g ) is the creation (re- 
spectively annihilation) operator of an electron in the 
effective medium a. We have therefore to determine 
two spin-dependent self-consistent effective baths, as in 
Ref . [2(J, but here we treat the full double exchange Hamil- 
tonian. Moreover, in the present study we go beyond the 
study of the paramagnetic state, since we are interested 
in the stability of the spontaneous ferromagnetic phase. 
The effective medium is characterized by the effective 
dynamical hybridization : 



J aa (oj) = 2J 



Kk\ 2 



i0+ 



£ f 

a,k.cr 



(8) 



The equation of motion for the Hamiltonian H cS gives 
the following equation for the retarded Green's function : 



G aa {u) 1 = cj - S aCT (o;) - J a a{u) 



(9) 



Comparing this equation to the following property that 
a Green's function satisfies on a Bethe lattice : 

G^iu)- 1 =u)-Z aa (u>)-t 2 G aa (cj) , (fO) 

we get the self-consistent equations of the DMFT : 

J acr (oj) =t 2 G aa {u) . (ff) 

As an impurity solver for the DMFT, we used the non- 
crossing approximation (NCA) 21 i 22 . Despite the known 
pathology of the NCA at very low temperature when ap- 
plied to the SIAM, it gives reliable results for tempera- 
tures down to a fraction of the Kondo temperature. As 
a self-consistent and conserving approximation, it also 
displays the correct scaling behavior and reproduces the 
relevant energy scales. In the DMFT framework the NCA 
has been applied successfully to various compounds like 




FIG. 1. Magnetization as a function of filling n for differ- 
ent temperatures (see legend), here J = 1.25. We remind 
the reader that the half- bandwidth is the energy unit. The 
ferromagnetic phase is centered around n = 1.4, at low tem- 
peratures it approaches n — 1. 



Lai-^Sr^TiO^and SrRuO^ allowing the 
first prediction concerning the orbitally selective Mott 
transition. 

Within the NCA, propagators P m (uj) and self-energies 
E m (w) for the 16 local eigenstates of the local Hamilto- 
nian H% are introduced. The sixteen coupled integral 
equations between local propagators and self-energies are 
solved numerically for each effective medium until con- 
vergence. 



IV. SPECTRAL PROPERTIES 

The phase diagram of this model was studied quite 
extensively by Peters et a^j 15 i 16 The phase diagrams ob- 
tained by our method are in good agreement with their 
results. The magnetization diagram in Fig. [T] is typi- 
cal for all magnetization diagrams for 0.75 < J < 2.0. 
At J = 0.5 there is no ferromagnetic phase at all, and 
at J = 0.75 its extension is small and centered around 
n ~ 1.4. The domain of ferromagnetism increases with J 
and for J > 1.5 the ferromagnetic phase includes quarter 
filling^ 

Since we did not introduce two sublattices, our numer- 
ical NCA procedure is unable to distinguish between or- 
bitally ordered and orbital liquid phases in the ferromag- 
netic domain. However, as will be explained in Sec. [Vj 
there are strong reasons to expect a saturated ferromag- 
netic phase with alternating orbital order near n = 1 at 
low temperatures (see also Ref. Hsl ) . For n = 1 + x elec- 
tron doping will destroy the orbital order (see Sec. IVI[) 
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FIG. 2. Spectral function for n — 1 — x and n — 2 — x, at 
temperature T = 4.3 x 10~ 3 and J = 1.0. The origin of the 
energy axis corresponds to the Fermi energy. 



in close analogy to the destruction of antifcrromagnctic 
long range order by hole doping in cupratc compounds. 
By the same analogy we expect a very small critical value 
of electron doping of only a few percent. Summarizing, 
we expect three different phases at zero temperature for 
n varying between 1 and 2 (and taking J = 1.25 as a 
representative value as in Fig.[T]): (i) The ferromagnetic 
and orbitally ordered phase close to n = 1, (ii) a sat- 
urated ferromagnetic and orbital liquid phase at higher 
values of electron doping, and (Hi) for n > 1.7 one gets 
a paramagnetic phase without orbital order. 

Figure [2] shows the spectral density p(e) = 
— ilm{G(e + iO)}. For n = 1 ± x or n = 2 ± x we can 
compare our results to those obtained in a recent paper 
by W.-C. Lee and P. W. Philips.™ Around n = 1 one 
clearly finds the series of three consecutive upper Hub- 
bard bands, with the spectral weight distribution 3:2:1 
(see Fig. [2]). We will see in Sec. IIV Bl that this spectral 
weight repartition is due to the double hopping term. 
Without explicit consideration of this term the three- 
peak structure of the unoccupied states is absent, lead- 
ing to qualitatively different behavior. Close to n = 2 
only two bands remain, each carrying the same spectral 
weight. The spectral weight transfer we observe between 
n = 1 — x and n = 2 + x is also in good agreement with 
earlier published analytic resultsJ^ The narrow peaks at 
the Fermi energy we observe in Fig. [2] are both quasi- 
particle peaks, we will discuss this in Sec. IIV Cl 




FIG. 3. Spectral function of the two spin states (red and blue) 
for different values of Hund's exchange J for n = 1.10, T = 
4.3 x 10~ :! and J = 0.5 - 1.25. (a) For J = 0.5 no magnetic 
phase is observed, (b)-(c) For J = 0.75 and J = 1.0 the 
magnetic phase does not include n = 1.10. (d) For J = 1.25 
the ferromagnetic phase is close to n = 1, here M = 0.68. 
The quasi-particle peak (inset) is discussed in Sec. IIV CI 



and the first upper band is U — As J increases the 

tendency towards ferromagnetism increases and is real- 
ized in Fig. [21(d) for J = 1.25. In the ferromagnetic case, 
the weight distribution 3:2:1 remains valid if we consider 
the sum over both spin directions, but its separation into 
up and down spin contributions is nontrivial. Because 
the system is slightly above n = 1 filling (Fig. © we also 
observe a spectral weight transfer to a band higher up in 
the spectrum, around e = 10, that will be filled only for 
n > 2. The insets show the quasi-particle peak that we 
will discuss in Sec. IIV CI 



B. Double hopping term 



A. Dependence on J 

It is instructive to study the J dependence of the spec- 
tral function (Fig.©. Near n = 1 (n = 1.10 for Fig.© 
the unoccupied part of the spectral function of the two- 
orbital Hubbard model is very characteristic with a spe- 
cific equidistant three peak structure. The energy differ- 
ence between the peaks is approximately 2 J which leads 
to a better separation of the peaks with the increase of J. 
The energy difference between the lower Hubbard band 



The pair hopping term is often neglected in studies 
investigating spectral functions at or around filling n = 1. 
This term is responsible for the splitting of a degenerate 
energy level (with energy U) at n = 2 into two different 
energy states U — J' and U + J' . These energy states are 
rather high in the energy spectrum and for < n < 1 
these bands are empty. In Fig. H]n = 1 + x, the Fermi 
energy lies on the edge of the first band of the series. 

The upper two graphs of Fig. [¥] characterize non- 
magnetic states of the system. The double hopping term 
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FIG. 4. J' = J (left) and J' = (right) for J = 1.0 (top) 
and J = 1.25 (bottom). When the double hopping term is 
included one finds back the spectral weight as calculated by 
Lee and Philips, when it is not included the spectral weight 
of the upper Hubbard bands changes. 



has influence only on the empty upper Hubbard bands, 
it does not change the magnetic phase of the system. In 
a system with the double hopping term we clearly find 
back the analytical results of Lee and Philips. Without 
this term the energy spectrum and the spectral weight 
in the upper spectrum change. It is the double hopping 
term that gives the spectral function its characteristic 
spectrum and spectral weight. 

In the ferromagnetic case (lower part of Fig. HJ, we 
can see that the spectrum and spectral weight change in 
much the same way as in the paramagnetic case. The 
series of three consecutive upper Hubbard bands with 
the spectral weight distribution 3:2:1 can be obtained by 
summing the spin up and spin down spectral weights. 
The magnetization changes slightly when one leaves out 
the double hopping term, but overall the phase diagram 
remains the same with or without the double hopping 
term. 



C. Quasi- particle peak 

For doping levels around n = 1 a quasi-particle peak 
is observed around the Fermi level (see inset of Fig. [3] 
and Fig. [5]). As our calculation gives fc-integrated Green 
functions, the small width of the peak Wqp indicates a 




FIG. 5. Quasi particle peak for n = 1.11 and n = 1.19, with 
J = 1.25 and T = 4.3 x 1(T 3 . The quasi-particle peak has the 
same spin as the majority spins of the lower Hubbard band. 

high effective mass m e ff ~ 1/Wqp of the corresponding 
quasi-particle. The peak is observed in the paramagnetic 
as well as in the ferromagnetic case. On the right hand 
side of Fig.[H]the width of the quasi-particle peak is larger 
than on the left hand side. That is due to two effects: the 
increase of magnetization and of doping. In general, the 
peak width increases if we pass from the paramagnetic 
to the ferromagnetic state. 

The quasi-particle peak is strongly doping dependent. 
It exists around quarter filling (n = 1) for electron and 
hole doping. The case of hole doping is shown in Fig. [2] 
(a). Exactly at n — 1, we have an insulating state and 
there is no quasi-particle peak which is presented in Fig. [5] 
(a). We compare different doping values (no doping, 20 
and 40 percent electron doping) and a tiny quasiparticle 
peaks on top of the main band exists only for low doping 
values. The high effective mass is due to short range or- 
bital correlations close to n = 1 as will be explained more 
in detail in the next chapters. For n = 1.40, the orbital 
correlations are weakened, and the quasi-particle peak 
merges with the main band. Similar spectral function 
had already been calculated before*^ It is interesting to 
note that a tiny quasi-particle peak exists also close to 
half-filling (see Fig. [U(b)). However, since we expect no 
tendency to orbital order at n = 2, that quasi-particle is 
probably of different origin and we did not investigate it 
further. 

In the ferromagnetic phase close to n = 1 the quasi- 
particle is of the same spin as the majority spin of the 
lower Hubbard band. That is in strong contrast to the 
antifcrromagnetic correlation between the local spin and 
the conduction electrons in the Kondo effect. Therefore, 
we have to find an alternative explanation that will be 
presented in Section IVT1 

V. EFFECTIVE HAMILTONIAN FOR n = 1 

To clarify the possible mechanism of quasi-particle for- 
mation one first has to analyze the possible phases at 
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FIG. 6. Appearance of the quasi particle peak for J = 1.25 
and different doping values. The magnetization is very close 
to one in all three cases. The tiny quasi-particle peak is only 
visible for n = 1.20, it is absent for n = 1, and it merges with 
the main band for high doping (n = 1.40). 
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FIG. 7. The four possible phases of the system at n = 1 are 
combinations of magnetic and orbital order. 



n=l. There are in total four possible phases of the 
system, fcrro- or antifcrromagnctic and alternating or 
homogeneous orbital order, see Fig [7] The ferromag- 
netic and homogeneous orbital case is highest in energy 
since hopping between neighboring sites is not allowed. 
For this reason this configuration can be excluded. The 
phase could also be antiferromagnetic with homogeneous 
or alternating orbital order. In this case hopping is pos- 
sible, but in the intermediate states the electrons meet 
at the same site with opposite spin and cannot benefit 
from Hund's exchange term. Lowest in energy is the fer- 
romagnetic state with alternating orbital order. 

To be more precise we derive now an effective low- 
energy spin-pseudospin Hamiltonian #le at n = 1. It 
corresponds to the low-energy limit of the two-orbital 
Hubbard Hamiltonian (P)-© assuming 



t<U,U', J <U and 



(12) 



Working in analogy to Refs. [l8|, [2fj, [27J, we obtain the 
effective Hamiltonian in terms of the spin S and orbital 



pseudospin f operators (see Appendix lAl for details): 



H 



R,9 



1 



U + J\4 
4 fl 



R R+g ' T R T R+g T R T R+g 



U-J\4 



r R T R+g 



(13) 



where we put U' — U — 2 J and J' = J to simplify the 
notation. This Hamiltonian corresponds to a special case 
a = 1 (a is the ratio between the transverse and longi- 
tudinal hopping terms) of the Hamiltonian of superex- 
change in the ab plane of alkali hyperoxide compounds 
R0 2 (R=K, Kb, Cs), cf. Eqs.(l)-(3) of Ref.13. 

It is not the aim of the present paper to analyze the 
low-energy Hamiltonian in detail. As it became clear 
from the qualitative considerations above, the ground 
state at n = 1 is the saturated ferromagnetic state with 
alternating orbital order. That is confirmed by the low- 
energy Hamiltonian. In the case of a saturated ferro- 
magnetic system we have Sr ■ Sr+ 9 = \. Then the low 
energy Hamiltonian (|13|) reduces to the form 



H 



2t 2 



LE 



U-3J 



E 



TR ■ TR + g 



R,g 

= JO TR ■ TR+g 

R,g 



(14) 
(15) 



that coincides with the Hcisenberg Hamiltonian, but in 
orbital space. Because Jo = 2< 2 /(C7 — 3J) is positive, 
the minimal energy configuration will be the one with 
alternating orbital order. The system is thus in a ferro- 
magnetic and alternating orbital phase. 



VI. ORBITAL POLARON 

The here derived low-energy Hamiltonian also gives in- 
formation about the mechanism that will be employed to 
make a charge carrier mobile on the lattice. The charge 
carrier can be a hole or an electron and we illustrate in 
Fig. [H] the case of electron doping (hole doping is equiv- 
alent). To describe the quasi-particle we have to add 
the kinetic energy term to the orbital Hcisenberg Hamil- 
tonian (|15p . In zeroth order approximation, this is the 
term ((IJ with charge fluctuations projected out (cf. the 
discussion after Eq. (|A7[) in Appendix |Aj The resulting 
model is the very well known t-J Hamiltonian. It was 
extensively studied in connection with high temperature 
superconductivity and it is quite interesting that our orig- 
inal model CE])-© of double exchange ferromagnetism can 
be reduced to the same form. The only difference is that 
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FIG. 8. Without orbital fluctuations the energy increases 
linearly with the length of the path traveled by the elec- 
tron, whereas with orbital fluctuations the electron can move 
around without increasing the energy proportionally to the 
length of the path. 



circular DOS corresponds to a Bethe lattice instead of a 
2D quadratic one. Also, it is not sure whether the single 
site DMFT is able to reproduce the correct quasi-particle 
band width. Nevertheless, it can be expected that the 
quasi-particle band width should be proportional to Jo 
close to n = 1. That is confirmed by the half- width at 
half- maximum of the quasi particle peak Wqp which can 
be read out of Fig. (JHJ) to be Wqp « 0.12 in a state 
which is close to saturated FM but still has low dop- 
ing. Therefore, Wqp is approximately equal to Jo which 
seems to be a realistic value. But one probably needs a 
more detailed method such as cluster DMFT to observe 
the correct dependencies of the quasi-particle bandwidth 
on the parameter values of the model for a given lattice 
type. 



VII. DISCUSSION AND CONCLUSION 



we have to replace spin fluctuations by orbital fluctua- 
tions. In the case of spin fluctuations, it is known that 
the one-hole state in the antiferromagnetic background of 
the t-J model gives rise to a coherent quasi-particle that 
is commonly called the spin polaron. Our case is very 
similar and leads to the orbital polaron. The orbital 
polaron was first proposed in the field of manganitcs. 
Taking a modified t-J model as the starting point, this 
kind of quasi-particle was shown to occur in the orbital- 
ordered ferromagnetic planes of LaMnOai^i Since then 
it was studied in several other worksii^— Here, we re- 
discover the orbital polaron in the complete double ex- 
change Hamiltonian including the double hopping term. 

We are not going to repeat the studies of spectral prop- 
erties in the framework of t-J Hamiltonians. We just re- 
mind here the basic arguments of the orbital polaron con- 
struction. The electron does not move by simply jumping 
from nearest neighbor to nearest neighbor since it would 
destroy the orbital order in such way (Fig. [H]). To be 
precise we rewrite Eq. (|15J) in a form that makes clear 
what the action of this low energy Hamiltonian is on the 
orbital pscudospin: 

#le = ffo„ + H Q± (16) 
= Jo £ f R f R+g + ^Y1 {j^R +a + ^i +a ) ■ 

R,g R,g 

It is the second term which allows for orbital fluctuations. 
Without it, the electron would be localized, since the 
application of H t leads to a string of misplaced orbitals 
with an energy increase that is proportional to the length 
of the path. It is the orbital fluctuation term Hqj_ which 
may reduce the length of the path. In such a way, it 
delocalizes the electron and we therefore expect a quasi- 
particle width of the order of Jo- 

The band-width of the spin polaron is known for the 
2D t-J model to be approximately 2 Jh where Jh is the 
exchange constant of the Heisenberg term. Our study 
is different in many respects. For instance, the semi- 



We investigated the complete two-orbital Hubbard 
model including the double-hopping term by means of 
the single site dynamical mean field theory (DMFT). We 
found ferromagnctism in the electron doped case n > 1 
that may extend up to n = 1 for sufficiently high values of 
the Hund exchange coupling J. The spectral function has 
two very characteristic features, a narrow quasi-particle 
peak at the Fermi level with a spin parallel to the major- 
ity spin direction and a characteristic equidistant three- 
peak structure with a weight distribution 3:2:1. We found 
the double-hopping term to be of crucial importance for 
obtaining the correct spectral function. 

Our study unifies several previous works. The mag- 
netic phase diagram of the present model, but without 
double hopping term, was already published bcforei 15 ! 16 
Our numerical results of the characteristic three-peak 
structure confirm the analytical predictions of Lee and 
Philips for the complete model including the double hop- 
ping termini We interpret the narrow quasi-particle peak 
as an orbital polaron that was known before in the field of 
manganitcs jii These former studies considered the one- 
hole problem in the effective, low-energy, t-J like Hamil- 
tonian. The present study is less restrictive since it treats 
the full two-orbital Hubbard model that is valid close to 
the Fermi level as well as for the high-energy features. 
Furthermore, it works for arbitrary doping. 

We admit a weakness of our approach since we did 
not calculate the two-sublattice problem. That would 
be necessary to obtain the correct alternating orbital- 
ordered phase at and close to n = 1 and for large U . It 
was obtained by Pruschke and Peters^ by DMFT and 
we confirm it here by the derivation of an effective spin- 
pseudospin Hamiltonian for n — 1. We have shown that 
this Hamiltonian leads to a ferromagnetic ground state 
with alternating orbital order. Adding additional charge 
carriers, we end up with the same t-J Hamiltonian as has 
been known for cuprate superconductors for a long time. 
This way we can interpret the observed quasi-particle as 
an orbital polaron in close analogy to the spin polaron. 
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The charge carriers have the tendency to destroy the or- 
bital order (see Fig. and we expect the orbitally or- 
dered phase in a small region around n = 1 only. That 
is the reason why we restricted our study to the orbital 
liquid phase where, however, short range alternating or- 
bital correlations persist. These correlations are visible 
in the spectral function by the small bandwidth of the 
quasi-particle. We have shown that the quasi-particlc 
peak broadens with increased doping indicating in such 
a way the decrease of orbital correlations. 

The two-orbital Hubbard Hamiltonian is relevant for 
several material classes. First of all it concerns man- 
ganitcs like LaMn03 and related compounds with an 



electron configuration close to t^g e g . In that case, 



(completely filled t2 g shells)^ — In undoped 
g orbitals are degenerate and there is no or- 



besides a partially filled e g shell that is well described 
by our model, an additional local spin (S = 3/2) ex- 
ists. That additional spin may lead to modifications that 
are outside the scope of the present model. It is bet- 
ter suited for nickelates like doped or undoped LiNiC>2, 
NaNi02 or AgNiC>2, with an electron configuration close 

to t% . 

LiNi0 2 the e 

bital or magnetic order in that specific two-dimensional 
triangular lattice (see Refs. [H and and references 
therein). However, electron doping leads to a ferromag- 
netic phase as it was observed recently^ and it reminds 
of our phase diagram (Fig.[T]). Manganites and nickelates 
are just two examples with partially filled e g shells but 
there exist many more possibilities in the field of tran- 
sition metal compounds (like LaVC>3 for example)^ As 
it was already mentioned in the introduction, the double 
exchange mechanism was also discussed for diluted mag- 
netic semiconductors (DMS)£ In many cases, the double 
exchange mechanism is in competition with other mech- 
anisms for ferromagnetism as for example itinerant mag- 
netism, direct exchange, Zener's kinetic p-d exchange, 
and so on. Our results for the spectral function provide 
a unique means to distinguish it from other mechanisms. 
This concerns especially the quasi-particle peak at the 
Fermi level. If, in a given ferromagnetic material a heavy 
quasi-particle peak with a spin parallel to the majority 
spin direction is observed, it would give a strong hint in 
favor of the double exchange mechanism. To observe it, 
we propose spin resolved photocmission with polarized 
photons^! Another characteristic feature is the three- 
peak structure with weight distribution 3:2:1 in the un- 
occupied part of the spectral function. Since the two 
distances between the three peaks are found to be equal, 
and very close to 2 J, it gives an experimental means to 
determine directly the Hund exchange coupling J. 

To interpret the quasi-particle peak, we concentrated 
the discussion on orbital fluctuations in a saturated 
ferromagnetic state. But our DMFT simulation is also 
able to describe spin fluctuations. These are present 
in the paramagnetic case. One might expect that the 
coupling to spin fluctuations makes the quasi-particle 
even more heavy. And indeed, we observed a reduction 
of the quasi-particlc band width if the system switches 



from the ferromagnetic to the paramagnetic phase. 
But a detailed theory of that process has still to be 
developed. The spin-pscudospin Hamiltonian which was 
derived in the present work may serve as a starting point 
for such a theory. 



Appendix A: Effective Hamiltonian 

Here we give the details of the derivation of the effec- 
tive low-energy Hamiltonian (TIB"]) , We write 



H d = J2 I*) Ei 



= oi, ( A1 ) 



where the eigenvalues Ei and eigenvectors \i) of the un- 
perturbed Hamiltonian H d are assumed to be known. 
Without loss of generality, we may assume that the per- 
turbation H t contains only non-diagonal terms. The op- 
erator 



W 



i¥=3 



Ei — Ej 



01, 



has the property 



H d ,W 



H,W - WH d 



-Ht 



(A2) 



(A3) 



Then, up to second order a canonical transformation 
gives an effective Hamiltonian 



H eS = cxp{-W)H exp(IF) = H d + H LE 
1 



(A4) 



Hi 



H,,W 



An explicit calculation gives 

#LE = \ E \f) tfjtjiD 



fji 



with 



D fji 
Ei4 



E 



fj 



E 3 . 



(A5) 

(A6) 
(A7) 



It is clear from Eqs. (|A5|) - (|A7[) that the transforma- 
tion may eliminate via Eq. (|A3[) only those terms in 
H t that connect the states with different energy, and 
it makes sense only when <C E^j. For our system 
this means that we eliminate only hoppings between the 
states which differ by total numbers of n-occupied sites, 
n = 0,1,2,3,4. At quarter filling the states |i) and |/) 
in Eq. (|A5[) contain one particle at each site, and in the 
intermediate states \ j) one site is empty and a neighbor- 
ing site is doubly occupied. The Hamiltonian i?LE in 
Eq. (|A5[) will contain the interactions between neighbor- 
ing sites. 
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In the Table I of Ref. [l7| a complete list of the cigen- 
states and eigenenergies of the Hamiltonian Hd ©-(HI) 
is given. We shall denote them as following: |0) is the 
vacuum state, number of electrons n = 0, energy Eq = 0. 
d\ a (J |0) is a one-particle state n = 1, with E a , a = 0; 
then, e.g., 



(A8) 



We will characterize two-particle states at lattice site r by 
I 



their energy and, when needed, the total spin projection 



r,U'- J,2a) = 4 iO X&,J >> 



r, U - J') = ± (4,a,t<«4 - <6,t<u) 1°) ' 

r > ^ + J ') = ~\ ««,t<-4 + <6,t<<a) 1°) • 



The effective second-order Hamiltonian (|A5|1 has the 
same form for every particular bond (a couple of neigh- 
boring sites) R, R' = R+g . So, we derive it for a two-site 
system, and the total Hamiltonian will be the sum over 
all bonds. The action of the operator H t on |i) 1IA8I) is 



S<j' -o [fia',ad Ra a d R,a,a + S a ' t - a d R d R 



(A9) 



Htd^ a J R ,^,. a , |0) = -t{5 a ,,J a ,,-J R _ a J R ^ a 

+ (R^R')\\0) 

= ij^cr'.cr^a'.-a (Sua - $ab) \R, U' - J, 2d) + 8 a ,- a -j= [8 a > , a (^<rf - (Saa + 5 a b) \R, U + J') 

+ 8a', a (6<rt - S^) {5 aa - 8 ab ) \R, U - J') + 6 a ' -a (Sat _ <^l) (Saa + S a b) \R, U' + J) 

+8*,-* (5 CTt + 8^) (5 aa - 5 ab ) \R, U' - J, 0)] + (R i?') }. 



I 

This expression provides us the matrix elements, which enter -&le (|A5|) . Now, we are ready to write the effective 

Hamiltonian for one bond and at n = 1. 

I 



#le = J2 d R,a,a d R> ,a> ,a> (°l d R >, a >, a >d R , a}(T H t \j) (j\ H t d R (j a d Ri p , s , |0) - ^ E d WJj/ . s ,d 



a, a' , {3,0' o,a' ,s,s f j 



l 



V — J 



l R,a,a u R'-a,a 



dR',-p, s d Rt p iS 5 as (8 aa — S a b) ($pa ~ Spb) 



7 d R,a,a d R',a,-a d R',P,s d R,P, S (Sat - <W ( S st ~ ^4-) 



U + J' 



+ jj Z j, d R, a ,a d R\a~a d R',l3-sdR,0,s (S aa ~ S ab ) (8p a - fy 6 ) (£4 - (5^) (£ st - S sl ) 



+ jj, , j d R, a> a d R',- a ,-a d R' ,sdR,0,s (Sat ~ S al) (<W ~ <W) 
1 

dR',-p,-sd Rt p j3 (5 aa — Sab) (Sfja ~ 8f3b) 



jjt y R,a,a R' , — a, — a 



(A10) 



Note that intermediate doubly occupied states \j) may be When n = 1, we may express the product of operators 
located on R as well as on R' sites. This gives a factor 2 at the same site d R a a d R p s , depending on the indices 
in the Ea.lfATOl). 
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a, P and a, s, via the product of spin S and pscudospin 
f operators, according to the rules (sec Eq. (3) of the 
Ref.Hl) 

{a, 13 = a) -> l/2 + af z , (cr, s = a) -> 1/2 + tr 5* 2 , 

(a,6)^f+, (t,i)^^ + , (All) 

(M)->f-, a,t)->£- 
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t 2 ^( 2 (1 

H ^---2^\jjr—jU 



R.g 

1 fl 



— 2?R ■ TR+g - + SrSr +3 



1- T R.TR+g + TrT r+ ^J 
J/- J' (2 + 2t R T R+9 ~ T R T R+g ~ T R T R+g S j 
fjl + j ( 2 _ 2l ~R T R+g + T R T R+g + T R T R+g 



Itlr 
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